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Average Structures of a Single Knotted Ring Polymer 
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Two types of average structures of a single knotted ring polymer are studied by Brownian 
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dynamics simulations. For a ring polymer with N segments, its structure is represented by a 
3iV-dimensional conformation vector consisting of the Cartesian coordinates of the segment 
positions relative to the center of mass of the ring polymer. The average structure is given 
by the average conformation vector, which is self-consistently defined as the average of the 
conformation vectors obtained from a simulation each of which is rotated to minimize its 
distance from the average conformation vector. From each conformation vector sampled in 
a simulation, 2N conformation vectors are generated by changing the numbering of the 
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segments. Among the 2N conformation vectors, the one closest to the average conformation 
vector is used for one type of the average structure. The other type of the averages structure 
uses all the conformation vectors generated from those sampled in a simulation. In the case of 
the former average structure, the knotted part of the average structure is delocalized for small 



N and becomes localized as N is increased. In the case of the latter average structure, the 
average structure changes from a double loop structure for small N to & single loop structure 
^vq , for large N, which indicates the localization-delocalization transition of the knotted part. 
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1. Introduction 



ulations, knot localization 



One of the important problems of the polymer physics is the effect of the topological con- 
straints on polymer systems. In the case of systems of linear polymers, the topological con- 



straints caused by the entanglement of the polymers vary temporally and the time-averaged 
properties of the constraints have been studied. 1 ' 2 In contrast, in the case of a single ring 
polymer, the topological constraints caused by its entanglement with itself, which is called a 
knot, are determined by the type of the knot and do not change with time. Therefore, a single 
ring polymer system can be considered as an ideal system for the study of topological effects 
and its investigation is expected to provide a basis for further understanding of the topological 
effects on polymer systems. The interest in single knotted ring polymers has been growing in 
recent years. 3-14 A crucial but unresolved problem in a single knotted ring polymer is how 
topological effects of the knot depend on the polymer length. There have been studies predict- 
ing that the topological effects vanish and properties of the knotted ring polymer become the 
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same as those of the unknotted ring polymer for infinitely long chains. 3-7 This prediction can 
be understood from the argument that the knotted part of a ring polymer becomes localized to 
a part of the polymer as the polymer length increases and that the rest of the polymer behave 
like an unknotted ring polymer. 7-14 The localization of the knotted part has been studied by 
measuring the chain length I of the knotted part through simulations. 8-12 It has been found 
that the equilibrium average (I) of I behaves as (I) oc iV* with t ~ 0.75 for a single knotted ring 
polymer of N segments in good solvent. 8 ' 9 The fact that < t < 1 indicates the localization of 
knotted part. It has been discussed and shown to be consistent with the simulation data that 
the behavior of (I) appears in the iV-dependence of the mean square radius of gyration as a 
leading correction to its scaling behavior. In contrast, a behavior (I) oc N, which corresponds 
to the knot derealization, has been found for a single knotted ring polymer in poor solvent. 
Similar results has been obtained for single ring polymers in two dimensions 12-14 and linear 
polymers. 10 Recently, the present authors studied the relaxation rate spectrum of single ring 
homopolymers in good solvent. 7 Because a ring homopolymer has the translational symmetry 
along the polymer, its relaxation rates are classified by a wave number q. In ref. 7, the relax- 
ation rate X q for a wave number q is estimated from the time t dependence of the correlation 
function C q (t) = iV -1 £\ £\ (1/3) ({n (t) - r c (t)) • (r 3 (0) - r c (0))} exp [-i27rg(t - j)/N] on 
the basis of the relaxation mode analysis, 15-17 where T{ (t) denotes the position of the ith seg- 
ment of a ring polymer with TV" segments at time t and r c (t) = ^2iL\ r i (t) is the position 
of the center of mass of the polymer. It is found for single ring polymers with the trefoil knot 
that the slowest relaxation rate for each N is given by X q with q = 2 for small values of N 
and that with q = 1 for large values of N. This transition is considered to correspond to the 
change of the structure of the ring polymer caused by the localization of the knotted part. In 
the studies mentioned above, however, the structures of single knotted ring polymers in three 
dimensions are not observed directly. 

The purpose of the present paper is to confirm the localization of the knotted part of a 
single knotted ring polymer by directly observing its average structure in three dimensions 
through simulations. The structure of a ring polymer with N segments is represented by a 
3iV-dimensional conformation vector consisting of the Cartesian coordinates of the segment 
positions relative to the center of mass of the ring polymer. The average structure is given 
by the average conformation vector, which is self-consistently defined as the average of the 
conformation vectors obtained from a simulation each of which is rotated to minimize its 
distance from the average conformation vector. 18 The average structure has been frequently 
used in simulation studies of biopolymers. Especially, static properties of fluctuations from 
the average structure have been studied by the principal component analysis method. 19-21 
Recently, the dynamic properties, such as the relaxation modes and rates, are also studied by 
the relaxation mode analysis method. 22 In the case of the studies of biopolymers, polymers 
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studied are usually heteropolymers. In contrast, a ring polymer studied in the present paper 
is a homopolymer and has the translational symmetry along the polymer chain. In the present 
paper, by utilizing the translational symmetry, we propose an extension of the definition of 
the average structure and examine how the localization of the knotted part appears in the 
new and conventional types of the average structures. 

The present paper is organized as follows. In §2, a model used in the present study and 
the two types of the average structures of single homopolymers are explained. The results of 
the simulations are presented in §3. Summary and discussion are given in the last section. 

2. Model and Average Structures 

In order to study a single knotted ring polymer in good solvent, Brownian dynamics 
simulations of a bead-spring model are performed. The dynamics of a single ring polymer 
with TV segments is described by the Langevin equation 

at C dr. 

Here, Tj(t) is a three-dimensional column vector consisting of the Cartesian coordinates of 
the position of the iih segment at time t and £ is the friction constant. The random force Wi 
acting on the ith segment is a Gaussian white stochastic process satisfying (wi jCe (t)) = and 
{wi,a{t)wjfi{t')) = (2/cbT/(") 6i,j5 ay p5{t — t'), where Wi >a , /«b and T denote the a-component of 
Wi, the Boltzmann constant and the temperature of the system, respectively. The potential 
F({rj}) = V{r\, . . . , rjy) represents the interaction between the segments. In the present 
paper, we use the potential given by 7 > 17 > 23 > 24 

N i-l N 

mrj}) = E E v ^ - r i\) + E ^(kw - rii), (2) 

i=2 j=l i=l 

where tn + i = n in the last summation of the right-hand side because the Nth segment is 
connected to the first segment. Here, Vr is given by the repulsive part of the Lennard-Jones 
potential 
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, , for r < 2 e a, 

V R (r) = <( LVr/ \rJ 4J ~ i (3) 

for r > 2e a, 

and represents the excluded volume interaction between all the segments. The potential Va, 
which is called a finitely extensible nonlinear elastic (FENE) potential, is given by 
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for r < Rn, 

(4) 

oo for r > Rq, 



and represents the attractive interaction between neighboring segments along the ring poly- 
mer. The same model is used for simulations of a single linear polymer with N segments, 
where the upper limit of the last summation of the right-hand side of eq. (2) is N — 1. 



3/15 



J. Phys. Soc. Jpn. Full Paper 

In the following, the average structures are explained for a linear or ring homopolymer 
with N segments. The average structures are estimated by using conformations of the polymer 
sampled in simulations. In a simulation, M conformations are sampled at interval of time AT 
after an initial equilibration time T\. Let Ri{m) denote the position of the iih segment relative 
to the center of mass of the polymer in the mth sample: 

Ri(m) = ri(t m ) - r c (t m ), (5) 

and t m = 7] + (m — I) AT, where m = 1,...,M. A structure of the polymer 
in the mth sample is represented by a 3iV-dimensional conformation vector C(m) = 
t ( t Ri(m), t R2(m), . . . ^Rpr^m)). From each C(m), K conformations, C(m,k) with k = 
1, . . . , K, each of which has the same statistical weight as C{m), can be generated by chang- 
ing the numbering of the segments, which is possible for homopolymers. Here, K = 2 for 
linear homopolymers and K = 2N for ring homopolymers. In both cases, C(m, k) is given by 
C(m, k) = t ( t i? n(lifc) (m), t i? n(2)fc) (m), . . . , t R n ^ Njk )('m)) > where n (h k) = (1 - 2r)ij + r(N + 1) 
with ij = ((i mod N) + 1, j = [(k - l)/2j , and r = (k — 1) mod 2 = k - 1 - 2j. 

Note that ij represents the integer between 1 and N which is equal to i + j modulo N. For 
odd or even values of k, n(i, k) is equal to ij or N — ij + 1, respectively. 

From the set of conformation vectors {C(m, k);m = 1, . . . , M and k = 1, . . . , K}, the 
average conformation vector C av * = t ( t i?f v * , t i?| v * , . . . , l RJ^*) , which gives one of the two types 
of the average structures, is calculated as follows. Here, the superscript av* denotes either 
avO or avl represeinting the type-0 or type-1 average structure explained in the following, 
respectively. For each conformation vector C{m, k), a three-dimensional rotation TZ m ,k which 
minimizes the square distance d 2 {m,k) between the conformation vector lZ m ,k (c(m,k)j 
obtained by the rotation of C(m, k) and the average conformation vector (7 av * is determined. 18 
Here, K mM (c(m, kfj = ^Tl m}k (R n(ljk) (m)) ^Tl^ (R n ( 2 ,k)(m)) , . . . ^Um^ {Rn{N,k){™))) 
and 

2 N 

d 2 (m, k) = (ll m , k (C{m, k)) - C av *) = Y, (K m ,k {Rn(i,k)(m)) - RT*)" ■ (6) 

i=i 

Thus, the rotation 1Z m ^ fits the conformation C{m,k) to the average conformation C av *. 

In the case of the conventional average structure, which has been used for heteropolymers, 
the ith segment of a sampled structure is fitted to the ith segment of the average structure. 
We call this type of average structure the type-0 average structure (type-0 AS). Because the 
numbering of the segments is crucial, the conformation vectors C(m,k), k = 1,...,K are 
considered to represent different structures, although they are obtained from C{m) by only 
changing the numbering of the segments. The average conformation vector (7 av0 for the type-0 
AS is then given by 

M K 
m=l k=l 
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In the following, we propose a new type of an average structure for homopolymers, which 
we call the type-1 average structure (type-1 AS). Because all the segments of a homopolymer 
are equivalent, the conformation vectors C(m, k), k = 1, . . . , K are considered to represent the 
same structure as C(m), although they have different numbering of the segments. Therefore 
the best fit of a sampled conformation C(m) to the average conformation is given by the 
rotation 'R- m) k miYl (m)i where A; m i n (m) is the value of k which gives the smallest value of d 2 (m, k) 
with m fixed. Thus, the average conformation vector (7 avl for the type-1 AS is given by 



Because the definition of the rotation lZ rn ^ contains C av *, eqs. (7) and (8) should be 
solved self-consistently. In practice, (7 av * is calculated iteratively. By using the nth candidate 
for the average conformation vector C^ v * calculated from the nth iteration, which may be 
chosen as one of C(m, k) for the first iteration (n = 0), all the rotations 1Z m ,k are determined. 
Then, the right-hand side of eq. (7) or (8) is calculated and the result is used as the next 
candidate C^Yi- The calculation is iterated for n = 0, 1, 2, . . . until the distance between C^* 
and C^+i becomes sufficiently small and the self-consistency of eq. (7) or (8) is achieved. 

3. Results 

Brownian dynamics simulations of the model described in the previous section are per- 
formed for a single linear polymer and single ring polymers with the trivial knot and the trefoil 
knot. The following parameters are used: 7 ' 17,23 k-^T/e = l,ka 2 /e = 30 and Rq/ct = 1.5. The 
Euler algorithm with a time step At = 10~ 4 Cc 2 /e is employed for a numerical integration of 
the equation of motion (1). Hereafter, we set a = 1, ( = 1 and e = 1. 

In order to calculate the average structures, conformations of a single polymer with N 
segments are taken from a simulation every AT = 10~ 2 t(N). For single linear polymers, the 
time t(N) is chosen as t(N) = A-^N XL , which corresponds to the behavior of the longest 
relaxation time of single linear polymers. 1,2 ' 16,17 For single ring polymers, t(N) is chosen as 
t(N) = ArN Xr , which corresponds to the behavior of the longest relaxation time for the 
wave number q = 1 of single ring polymers with the trivial knot. 7 The parameters are chosen 
as Al = 21.05, xl = 2.22, Ar = 22.96 and xr = 2.095. The values for linear polymers 
are estimated from the relaxation mode analysis of the time-displaced correlation matrix 
((fi(t) — r c (t)) ■ (rj(0) — r c (0))) calculated through simulations. The values for ring polymers 
are taken from the previous study. 7 The equilibration time Tj is chosen as T| ~ Wt(N). 
The number of sampled conformations M is given by M = 10 5 for N = 30, 40, 60 and 80, 
M = 2 x 10 5 for N = 110 and 120, M = 1.5 x 10 5 for N = 160 and M = 5 x 10 4 for N = 240, 
respectively. 

Figure 1 shows the average structures for a single linear polymer with N = 40 segments. 
In this and the following figures, the bonds connecting adjacent segments of a polymer are 





m=l 
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represented by cylinders and the center of mass of the polymer is shown as a sphere. The 
axes of the Cartesian coordinates are chosen to be the principal axes of the moment of in- 
ertia tensor and the length scale is chosen to normalize the contour length of each average 
structure. In the figures, only the directions of the axes are indicated. The origin of the 
Cartesian coordinates is chosen to be the center of mass of each average structure in the 
following descriptions. The type-0 AS shown in Fig. 1(a) has a parabolic shape in the x-y 
plane, which has C2 symmetry about the y axis. This result agrees with that of the previous 
study of the average structure of single linear polymers. 25 The C2 symmetry is induced by 
the symmetry of a single linear homopolymer that the ith and (N — i + l)th segments are 
equivalent, because the numbering of the segments is conserved in the type-0 AS. On the 
other hand, the type-1 AS shown in Fig. 1(b) does not have the symmetry, although it has 
a similar shape in the x-y plane. Each sampled conformation C(m) does not have the C2 
symmetry, that is, there is no rotation which changes C(m) = C(m,l) into C(m, 2). This 
asymmetry in the sampled conformations is conserved in the type-1 AS, because it neglects 
the numbering of the segments. The reason why the type-1 AS is in the x-y plane is as 
follows. Let A4 denote a reflection operation with respect to a plane containing the origin, 
which is the center of mass of conformations in the present case. For a rotation operation 
TZ, we define another rotation operation as 1Z' = A41ZA4. Then, (iZ i^C{m,k)^ — C av *^ = 

(M (n(c(m,k)^ -M{C av *)y = (n' (M(c(m,k)^ -M{C W *)^ holds. Therefore, 
if TZ m ,k min (rn) gives the smallest square distance calculated from C(m) and C avl , lZ' m k . ^ 
gives the smallest square distance calculated from M (C(m)) and M (C avl ). By considering 
the mirror image of eq. (8), we have 

1 M 

m=l 
1 M 

= M E {M (C(m, fc min (m)))) . (9) 

m=l 

This equation means that if the type-1 AS calculated from M samples C(m), m = 1, . . . ,M 
is given by C avl then that calculated from M samples M (C(m)), m = 1, . . . ,M is given 
by M (C avl ). For a sampled conformation C(m), its mirror image M. (C(m)) has the same 
statistical weight as C(m) in the ensemble of conformations, because a single linear polymer 
has no chirality. In other words, the ensemble of conformations M (C(m)) is the same as 
that of conformations C(m). Therefore, the average conformation C avl and its mirror image 
M (C avl ) represent the same conformation, that is, there exists a rotation 1Z which satis- 
fies C avl =TZ(M (C av1 )). It should be a conformation in a plane or a three-dimensional 
conformation with a plane of reflection symmetry. 

Figure 2 shows the average structures of a single ring polymer with the trivial knot for 
the case of iV = 40. The type-0 AS shown in Fig. 2(a) forms a regular polygon of iV = 40 
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Fig. 1. The type-0 AS (a) and the type-1 AS (b) of a single linear polymer with N = 40. In each 
figure, the cylinders and the sphere represent the bonds connecting adjacent segments and the 
center of mass of the polymer, respectively. 

sides in the x-y plane. As mentioned before, 2N conformations C(m, k), k = 1, . . . , 2N are 
generated from a sampled conformation C(m) by changing the numbering of the segments. 
The 2N ways of changing the numbering correspond to the 2N symmetry operations of the 
dihedral group Djv- Because all the 2N conformations are used with the same statistical 
weight in the calculation of the type-0 AS, the type-0 AS should be invariant under the 2N 
ways of changing the numbering. Therefore, the type-0 AS should be a regular polygon of 
N sides in two dimensions, which has the Djy symmetry, if all the position vectors in the 
average structure are different each other. In contrast, the type-1 AS shown in Fig. 2(b) forms 
a distorted polygon, which does not have the D;v symmetry, in the x-y plane. The form of 
the type-1 AS is explained in the same way as before. Because the numbering of the segments 
is neglected in the type-1 AS, the deviations from the D^r symmetric form in the sampled 
conformations are conserved in the type-1 AS. Because a sampled conformation and its mirror 
image have the same statistical weight in the ensemble of conformations of a ring polymer 
with the trivial knot, which has no chirality, the average structures should have reflection 
symmetry. A structure in a plane is one of the possible structures. 

Figure 3 shows the type-0 ASs of single ring polymers with the trefoil knot for N = 30, 
40, 60, 80, 110, 120, 160 and 240. As explained before, the type-0 AS should be invariant 
under the symmetry operations corresponding to the 2N ways of changing the numbering 
of the segments. Therefore, the deviations from the x-y plane in the average structures for 
80 < N < 160 are considered to be due to statistical errors. If the deviations for 80 < N < 160 
are ignored, the type-0 ASs shown in Fig. 3 form regular polygons in the x-y plane. In order 
to see how the segments are distributed on the regular polygon, we calculate the angle 6i 
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between the position vectors projected onto the x-y plane of the first and the ith segments: 

oavO oavO i oavO oavO 
COS#j = = (10) 



jjavO c>av0 pavO pavO 

sin ft = =^Sg (11) 

v^:*° 2 + R rf^Rt:f + ^° 2 

Here, -Rf^ 1 denotes the a-component of the position vector Rf v0 of the ith segment relative 
to the center of mass in the type-0 AS. Figure 4 shows a plot of Oijlix versus (i — 1)/N. It 
is clearly seen that the type-0 AS forms a double loop for A < 110 and a single loop for 
N > 120. In other words, the type-0 AS forms a regular polygon of N/2 sides for N < 110, 
where N is even, and a regular polygon of N sides for N > 120. Note that in the case of the 
double loop structure with even N, the symmetry operations corresponding to the 2A ways of 
changing the numbering of the segments become those of T) N / 2 , because Rf v0 = ^ or 
i = 1, • • ■ , A/2. The region of the double loop structure and that of the single loop structure 
must be separated by a transition segment number A t av0 between 110 and 120. This transition 
is considered to correspond to the localization-delocalization transition of the knotted part, 
because if the knotted part of a conformation 1Z m ,k (c(m,k) \ is delocalized, its projection 
onto the plane of the type-0 AS encompasses the center of mass while the projection of the 
knotted part does not encompass the center of mass if the knotted part is localized. Note that 
the form of the type-0 AS of the ring polymer with the trefoil knot for N > 120, which is a 
regular polygon of N sides, is the same as that with the trivial knot. 

Figure 5 shows the type-1 ASs of single ring polymers with the trefoil knot for N = 30, 
40, 60, 80, 110, 120, 160 and 240. In the case of the trefoil knot, the type-1 ASs do not 
have reflection symmetry. This is because a ring polymer with the trefoil knot has chirality 
and therefore the mirror image of a sampled conformation C(m) is topologically different 
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from C(m) and does not appear in the ensemble of conformations to which C(m) belongs. 
Moreover, the knot type of the type-1 AS is identical with that of the original structure for 
each N. It is clearly seen that the knotted part is delocalized for small values of N and is 
localized for large values of N. The crossover segment number iVJ vl ~ 120 from the delocalized 
state, where the center of mass is inside of the knotted part, to the localized state, where the 
center of mass is outside of the knotted part, is consistent with the transition segment number 
iV t av0 of the type-0 AS. 

4. Summary and Discussion 

In the present paper, the two types of average structures are calculated for a single linear 
polymer and single ring polymers by Brownian dynamics simulations. The average conforma- 
tion vector, which specifies the average structure, is self-consistently defined by eqs. (7) and 
(8): It is the average of the sampled conformation vectors each of which is rotated to minimize 
its distance from the average conformation vector. From a conformation vector C{m) of single 
homopolymers, K conformation vectors C(m, k) with k = 1, . . . , K are generated by changing 
the numbering of the segments, where K = 2 for linear homopolymers and K = 2N for ring 
homopolymers with N segments. In the calculation of the type-0 average structure, which 
has been used for heteropolymers, the position of the ith segment of a sampled conformation 
is fitted to that of the ith segment of the average conformation. Therefore, in the case of 
homopolymers, the above-mentioned K conformation vectors represent different structures 
and are used with the same statistical weight in the calculation as shown in eq. (7). Thus, the 
type-0 AS should be invariant under the K ways of changing the numbering of the segments. 
For a single linear polymer, the type-0 AS has a parabolic shape in the x-y plane with the C2 
symmetry about the y axis, where the two symmetry operations of the group C2 correspond 
to the two ways of changing the numbering. In the case of a single ring polymer of iV segments 
with the trivial knot, the type-0 AS is given by a regular polygon of iV sides with the D/v 
symmetry, where the 2N symmetry operations of the group D^r correspond to the 2N ways 
of changing the numbering. 

In the calculation of the type-1 average structure, which is proposed in the present paper, 
differences in the numbering of the segments are neglected and all the K conformation vectors 
represent the same structure. Therefore, among the K conformation vectors, the only one 
conformation C(m, k m - m (m)), which becomes the closest to the average conformation vector 
by the rotation 7£ mi fc min (m), is used in the calculation of the type-1 AS as shown in eq. (8). 
Thus, deviations from the symmetric form of the type-0 AS in the sampled conformations are 
conserved in the type-1 AS. The type-1 AS of a single linear polymer has a distorted parabolic 
form without the C2 symmetry in the x-y plane. For a ring polymer of N segments with the 
trivial knot, the type-1 AS forms a distorted polygon of iV sides without the D^r symmetry 
in the x-y plane. It is shown that the mirror image of the average structure calculated from 
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conformations of a polymer gives the average structure calculated from the mirror images of 
the conformations. In the ensemble of conformations of a linear polymer or a ring polymer 
with the trivial knot, which has no chirality, a conformation and its mirror image have the 
same statistical weight and the ensemble of the mirror images of the original conformations 
is the same as the original ensemble. Therefore the average structure and its mirror image 
represent the same structure. Thus, the average structure is given by a conformation in a 
plane or a three-dimensional conformation with a plane of reflection symmetry. All the average 
structures of a single linear polymer and single ring polymers with the trivial knot are planar 
and consistent with the above consideration. 

The type-0 AS of a single ring polymer of N segments with the trefoil knot forms a 
double loop on a regular polygon of N/2 sides for even N < 110 and a single loop on a 
regular polygon of N sides for N > 120. These structures are invariant under the symmetry 
operations corresponding to the 2N ways of changing the numbering. The transition from the 
double loop structure to the single loop structure occurs at the transition segment number 
iV t av0 between 110 and 120. This transition is considered to correspond to the localization- 
delocalization transition of the knotted part. Because a ring polymer with the trefoil knot 
is chiral, its type-1 AS does not have reflection symmetry and has the same knot type as 
the sampled conformations. The knotted part of the type-1 AS is expanded along the whole 
structure for small values of N and localized to a part of the structure for large values of N. 
The crossover from the delocalized state to the localized state occurs around the crossover 
segment number N£ vl ~ 120, which is consistent with the value of iV t av0 . The transition of 
the type-0 AS and the crossover of the type-1 AS furnish strong evidence for the localization 
of the knotted part predicted in the previous study, where the transition segment number is 
between 120 and 160. 7 

It is demonstrated that the analysis of the average structures is useful for studying the 
localization of the knotted part of single knotted ring polymers. The present paper provides 
strong evidence for the knot localization by directly observing the average structures in three 
dimensions. In the case of single ring polymers with the trefoil knot, the transition of the 
type-0 AS from the double loop structure to the single loop structure is observed as the 
number of the segments N is increased, which is considered to correspond to the localization- 
delocalization transition of the knotted part. The trefoil knot is a (3, 2)-torus knot. Here, a 
(p, g)-torus knot is given by a curve on the surface of a torus which winds p times around the 
center line of the torus and revolves q times along the center line of the torus. 26 In the case of 
the trefoil knot, the curve revolves twice along the center line of the torus. Therefore, it seems 
natural that the type-0 AS forms a double loop structure for small N and that the transition 
to a single loop structure occurs as N is increased. It is expected that the similar transitions 
occur in the type-0 ASs of single ring polymers with another torus knot, such as 5i, 7i, or 
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9i . It is interesting to study how the type-0 ASs of single ring polymers with a nontorus knot 
depends on N. The type-1 ASs of single ring polymers with the trefoil knot have the same 
knot type as the original conformations and show the knot localization directly. Single ring 
polymers with the trefoil knot have chirality. Therefore, their type-1 ASs need not to have 
reflection symmetry. Thus, single knotted ring polymers with chirality have less constraints 
on the shapes of their type-1 ASs. Because torus knots have chirality, it is expected that the 
type-1 ASs of single ring polymers with another torus knot show the knot localization in three 
dimensions. It is interesting to study what kind of structures are formed by the type-1 ASs 
of single knotted ring polymers without chirality. They must be planar structures or three- 
dimensional structures with reflection symmetry. The study of the average structures of single 
ring polymers with the figure eight knot is very important, because the figure eight knot is the 
next simplest prime knot and a nontorus knot without chirality. The study in this direction 
is in progress. 
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Fig. 3. The type-0 ASs of a single ring polymer with the trefoil knot for N = 30, 40, 60, 80, 110, 120, 
160 and 240. The meaning of the cylinders and the spheres is the same as in Fig. 1. 
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Fig. 4. A plot of 9i/{2ir) versus (i — 1)/N for the type-0 ASs of single ring polymers with the trefoil 
knot for N = 60, 110, 120 and 240. For clarity, only the data for (i - 1)/N = n/60 with n = 
0, 1, . . . , 59 arc plotted for N = 120 and 240. 
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Fig. 5. The type-1 ASs of a single ring polymer with the trefoil knot for N = 30, 40, 60, 80, 110, 120, 
160 and 240. The meaning of the cylinders and the spheres is the same as in Fig. 1. 
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